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ORIGINAL ARTICLE 



Multiple introductions of highly pathogenic avian 
influenza H5N1 viruses into Bangladesh 

Atanaska Marinova-Petkova\ Mohammed M Feeroz^, SM Rabiul Alam^, M Kamrul Hasan^, Sharmin Akhtar^, 
Lisa Jones-EngeP, David Walker^, Laura McClenaghan^, Adam Rubrum\ John Franks^, Patrick Seiler^, 
Trushar Jeevan\ Pamela McKenzie^, Scott Krauss^, Richard J Webby^ and Robert G Webster^ 

Highly pathogenic H5N1 and low pathogenic H9N2 influenza viruses are endemic to poultry markets in Bangladesh and have 
cocirculated since 2008. H9N2 influenza viruses circulated constantly in the poultry markets, whereas highly pathogenic H5N1 
viruses occurred sporadically, with peaks of activity in cooler months. Thirty highly pathogenic H5N1 influenza viruses isolated from 
poultry were characterized by antigenic, molecular, and phylogenetic analyses. Highly pathogenic H5N1 influenza viruses from clades 
2.2.2 and 2.3.2.1 were isolated from live bird markets only. Phylogenetic analysis of the 30 H5N1 isolates revealed multiple 
introductions of H5N1 influenza viruses in Bangladesh. There was no reassortment between the local H9N2 influenza viruses and 
H5N 1 genotype, despite their prolonged cocirculation. However, we detected two reassortant H5N 1 viruses, carrying the M gene from 
the Chinese H9N2 lineage, which briefly circulated in the Bangladesh poultry markets and then disappeared. On the other hand, 
interclade reassortment occurred within H5N1 lineages and played a role in the genesis of the currently dominant H5N1 viruses in 
Bangladesh. Few 'human-like' mutations in H5N1 may account for the limited number of human cases. Antigenically, clade 2.3.2.1 
H5N1 viruses in Bangladesh have evolved since their introduction and are currently mainly homogenous, and show evidence of recent 
antigenic drift. Although reassortants containing H9N2 genes were detected in live poultry markets in Bangladesh, these reassortants 
failed to supplant the dominant H5N1 lineage. 
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INTRODUCTION 

The highly pathogenic avian influenza (HPAI) A (H5N1) virus first 
appeared in 1996 in geese in Guangdong, China, and continues to 
circulate in the poultry population in many countries, mainly in 
Asia. Although the HPAI H5N1 virus has not yet demonstrated the 
ability to transmit efficiently between humans, it remains a public 
health concern because of the high case-fatality rates associated with 
H5N1 infection. Recent studies on ferrets show that limited changes 
are required for the H5N1 virus to acquire the ability for airborne 
transmission.^'^ 

Of the 32 distinguishable genetic clades of H5N1 viruses,^ those that 
are currently circulating are clades 1.1 (Cambodia and Vietnam); 
2.1.3.2 (Indonesia); 2.2.1 (Egypt); 2.3.2.1 (Bangladesh, China, India, 
Indonesia, Korea, Nepal and Vietnam); 2.3.4.2 (China); and 7.2 
(China and Vietnam).'' Clade 2.3.2.1 has been particularly associated 
with wild-bird infections,^' and after it was first detected in migratory 
birds in Mongolia in 2009,^ it spread to eastern Europe in 2010.^'^ 

Bangladesh has one of the highest densities of human (1072 people/ 
km^) and poultry populations (1194 birds/km^) in the world.** 
Approximately 75% of the human population in Bangladesh depends 
heavily on poultry as a source of meat and household income. Thus, 
there is a high risk of rapid disease spread and transmission to a large 



number of people in the event that a virus with a pandemic potential 
emerges in Bangladesh. 

The HPAI H5N1 virus was first detected in Bangladesh on 5 
February 2007, on a poultry farm in Jamalpur district, Dhaka division, 
and as of 1 November 2013, there have been 549 outbreaks reported in 
the country.^ From February 2007 until the end of 2010, the circulat- 
ing HPAI H5N1 viruses in Bangladesh were all from clade 2.2, sub- 
lineage III; however, at the beginning of 2011, new introductions of 
clades 2.3.2.1 and 2.3.4.2 were detected."*''' 

As of 1 November 2013, seven human cases of HPAI H5N1 infec- 
tion have been reported in Bangladesh.'^ All the cases had a history of 
exposure to poultry, including sick or dead birds. The first three cases 
occurred in 2008-201 1 in Kamalpur, Dhaka, in children under the age 
of two years, and the viruses isolated from the first and the second case 
belonged to clades 2.2 and 2.2.2, respectively. ''''^^ Three cases were 
reported in 2012 in adult males working at the live bird markets in 
Dhaka city and all three viruses belonged to clade 2.3.2.1.'' On 18 
February 2013, the seventh case and first human death caused by 
HPAI H5N1 was reported in a 23-month-old boy from Comilla. 
Influenza A (H5N1) virus from clade 2.3.2.1 was isolated from the 
patient's respiratory swab and from tissue samples from sick chickens 
with which the boy had been in contact.'^ 



^Department of Infectious Diseases, Sf. Jude Children's Research Hospital, Memphis, TN 38105, USA; ^Department of Zoology, Jahangirnagar University, Dhaka-1342, 
Bangladesh and ^National Primate Research Center, University of Washington, Seattle, WA 98195, USA 
Correspondence: RG Webster 
E-mail: Roberf.Webster@Sfjude.org 

Received 18 November 2013; revised 17 December 2013; accepted 18 December 2013 



H5N1 viruses in Bangladesli 

Marinova-Petkova etal 



In addition to the sporadic spread of H5N1 viruses, H9N2 influenza 
viruses have been isolated year-round from birds at the live bird mar- 
kets in Dhaka city since 2008.^^''^ Before these studies, isolation of 
H9N2 viruses in Bangladesh had not been reported. However, on the 
basis of molecular and phylogenetic analyses, Shanmuganatham 
etal}^ hypothesized that the H9N2 subtype had been introduced into 
Bangladesh from Pakistan through poultry movements in the early 
2000s. 

In this study, we performed antigenic, molecular and phylogenetic 
characterization of the dominant H5N1 lineages in Bangladesh from 
November 2009 to June 2013 and found high homology among the 
majority of H5N1 viruses. We also detected the separate introduction 
of two reassortant HPAI viruses of subtype H5N1 as well as interclade 
reassortment, but did not find reassortants with cocirculating H9N2 
viruses. 

MATERIALS AND METHODS 

Sample collection 

This study report results from ongoing active avian influenza surveil- 
lance in Bangladesh that started in November 2008. The target loca- 
tions for sample collection were five retail poultry markets and a pet 
bird market in Dhaka city; six poultry farms (chicken and duck farms) 
in the Dhaka and Rajsashi divisions; the Baikka Beel wetlands in 
Sreemangal, located in the Moulvibazar district (200 km northeast 
of Dhaka); and a lake in Savar, Dhaka, which serves as a feeding station 
for migratory wild birds. Samples from the retail markets and the 
farms were collected throughout the year, whereas environmental 



samples from the wetlands and the lake were collected only during 
the winter and spring seasons when the wild birds are present during 
migration. 

Sample screening and virus isolation 

All samples were screened for avian influenza viruses (AIVs) by real- 
time reverse transcription polymerase chain reaction, using influenza 
A-specific primers and probes^^ as previously described.'^ All influ- 
enza A real-time reverse transcription polymerase chain reaction-pos- 
itive samples were inoculated into 10-day-old embryonated chicken 
eggs to confirm the presence or absence of virus infection. 

Viral RNA was extracted from allantoic fluid, using the RNeasy 
Mini Kit (Qiagen, Valencia, CA, USA) and isolates were initially 
screened by sequencing for influenza viruses subtypes H5 and H9 
and the Newcastle disease virus, since all three viruses commonly 
circulate in Bangladesh. 

Antigenic analysis 

Preparation of post-infection ferret antisera for hemagglutination 
inhibition (HI) assays. All post-infection ferret antisera used in the 
HI assays in this study were produced in United States Department of 
Agriculture-approved biosafety level 3 enhanced facihties at St Jude 
Children's Research Hospital (St Jude, Memphis, TN, USA). Because 
the antisera a-RG-A/barn swallow/Hong Kong/1161/2010 (H5N1), 
a-A/environment/Bangladesh/15121/2012 (H5N1) and a-A/duck/ 
Bangladesh/19097/2013 (H5N1) induced very low antibody titers after 
infection, the immune response in ferrets was boosted with Freund's 



Table 1 Avian influenza A (H5N1) virus isolates from live bird markets in Bangladesh from December 2010 to March 2013 


Isolate 


Clade 


Isolation date (mm/dd/yy) 


Host species 


Sample type 


Retail market 


A/Ck/BD/9636/10 


2.2.2 


12/11/10 


Chicken 


C 


1 


A/Ck/BD/9675/11 


2.2.2 


01/18/11 


Chicken 


F 


2 


A/Ck/BD/12197/11 


2.3.2.1 


08/14/11 


Chicken 


W 


1 


A/Ck/BD/15078/12 


2.3.2.1 


01/18/12 


Chicken 


C 


3 


A/Ck/BD/15079/12 


2.3.2.1 


01/18/12 


Chicken 


OP 


3 


A/Ck/BD/15083/12 


2.3.2.1 


01/18/12 


Chicken 


OP 


3 


A/Ck/BD/15085/12 


2.3.2.1 


01/18/12 


Chicken 


OP 


3 


A/Env/BD/15105/12 


2.3.2.1 


01/18/12 


Quail 


F 


4 


A/Env/BD/151 14/12 


2.3.2.1 


01/18/12 


Quail 


F 


4 


A/Env/BD/15121/12 


2.3.2.1 


01/18/12 


Quail 


F 


4 


A/Env/BD/15126/12 


2.3.2.1 


01/18/12 


Quail 


F 


4 


A/Env/BD/15145/12 


2.3.2.1 


01/18/12 


Quail 


F 


4 


A/Env/BD/15147/12 


2.3.2.1 


01/18/12 


Quail 


F 


4 


A/Env/BD/15165/12 


2.3.2.1 


01/18/12 


Quail 


F 


4 


A/Ck/BD/15210/12 


2.3.2.1 


01/18/12 


Chicken 


OP 


1 


A/Env/BD/15661/12 


2.3.2.1 


02/08/12 


Chicken 


F 


2 


A/Env/BD/16298/12 


2.3.2.1 


03/31/12 


Chicken 


W 


2 


A/Dk/BD/17753/12 


2.3.2.1 


07/24/12 


Duck 


OP 


5 


A/Ck/BD/18061/12 


2.3.2.1 


09/07/12 


Chicken 


c 


5 


A/Ck/BD/18065/12 


2.3.2.1 


09/07/12 


Chicken 


c 


5 


A/Ck/BD/18248/12 


2.3.2.1 


10/21/12 


Chicken 


c 


5 


A/Env/BD/18948/13 


2.3.2.1 


01/28/13 


Duck 


w 


1 


A/Dk/BD/18949/13 


2.3.2.1 


01/28/13 


Duck 


w 


1 


A/Ck/BD/18984/13 


2.3.2.1 


01/28/13 


Chicken 


w 


2 


A/Env/BD/18989/13 


2.3.2.1 


01/28/13 


Chicken 


w 


2 


A/Dk/BD/19097/13 


2.3.2.1 


02/14/13 


Duck 


OP 


5 


A/Dk/BD/19099/13 


2.3.2.1 


02/14/13 


Duck 


OP 


5 


A/Quail/BD/19250/13 


2.3.2.1 


03/30/13 


Quail 


OP 


4 


A/Quail/BD/19254/13 


2.3.2.1 


03/30/13 


Quail 


OP 


4 


A/Env/BD/19338/13 


2.3.2.1 


03/30/13 


Chicken 


w 


1 



Abbreviations: BD, Bangladesh; C, cloacal; Ck, chicken; Dk, duck; Env, environment; F, fecal; OP, oropharyngeal; W, water. 

The five anonymous retail poultry markets In Dhaka city where the H5N1 viruses were detected, were assigned numbers 1, 2, 3, 4, and 5 respectively. 
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incomplete adjuvant (Invivogen, Inc., San Diego, CA, USA) at day 14 
post-infection. Initially, each of three groups of three ferrets was 
infected intranasally with one of the following hve viruses: 0.5 mL 10^ 
50% egg infectious dose (EID50) of wild-type A/duck/Bangladesh/ 
19097/2013(H5N1), 0.5 mL 10^ EID50 of wild-type A/environment/ 
Bangladesh/15121/2012 (H5N1) or 1.0 mL 10** EID50 of reverse genetics 
2-1-6 A/barn swallow/Hong Kong/1 161/2010 (H5N1), A/Puerto Rico/8/ 
34 (internal genes). Fourteen days post-infection, the surviving ferrets 
were bled and then boosted intramuscularly with 0.5 mL 10^ EID50 of 
the homologous virus, with Freund's incomplete adjuvant according to 
the instructions of the manufacturer. One ferret from the group receiv- 
ing A/environment/Bangladesh/15121/2012 (H5N1) did not survive the 
initial infection. 

HI assays. HI assays were done as previously described. The panel of 
antisera used in the HI assays included representatives from the cur- 
rently circulating genetic lineages of clade 2.3.2.1 in Asia and 
Bangladesh, as well as a representative from clade 2.2. Horse red blood 
cells (1%) were used for the assays. 

Sequencing and phylogenetic analysis. DNA sequencing was per- 
formed by the HartweU Center for Bioinformatics and Biotech- 
nology at St Jude, using universal primers''^ and subtype-specific 
primers (primer sequences available upon request), and the BigDye 



Terminator v3.1 Cycle Sequencing Kit (Invitrogen, Carlsbad, CA, 
USA) on 3730XL DNA analyzers (Applied Biosystems, Foster City, 
CA, USA), according to the manufacturer's recommendations. 

The hemagglutinin (HA) and neuraminidase (NA) genes of all 30 
H5N1 isolates from Bangladesh were sequenced, and 23 isolates from 
2010 to 2013 were selected for complete genome sequencing (Table 1). 
FuU-genome phylogenetic analysis of the sequenced Bangladeshi 
H5N1 viruses was performed, and representative nucleotide (nt) 
sequences of H5N1 clades that have been recently circulating in Asia 
and Egypt were retrieved from the National Center for Biotechnology 
Information (NCBI) Influenza Virus Sequence Database^" and the 
EpiFlu database of the Global Initiative on Sharing AH Influenza 
Data^'. All nucleotide sequences were aligned by using the ClustalW 
tool in MEGA 5 and trimmed to equal lengths, but at least 90% of the 
coding region as follows: HA, 1630 nt; NA, 1350 nt; polymerase basic 
protein 1 (PBl), 2213 nt; polymerase basic protein 2(PB2), 2214 nt; 
polymerase acidic (PA), 2069 nt; nucleoprotein (NP), 1416 nt; Matrix 
(M), 744 nt; and non-structural (NS), 693 nt. Phylogenetic relation- 
ships were inferred by the neighbor-joining method from 1000 boot- 
strap values, the topology was confirmed by the maximum likelihood 
method, and evolutionary analyses were conducted in MEGA 5.^^ 

The sequences obtained in this study were deposited in the Influenza 
Research Database and are available under GenBank accession numbers: 
AGZ62382, AGZ62383, KF874284-KF874289 and KF888395-KF888584. 




Figure 1 Samples collected and avian influenza viruses isolated in Bangladesh (November 2008-June 2013). Total number of samples collected was 19897 (3952 
oropharyngeal, 4113cloacal and 11 832 environmental samples, the last set consisting of 1044 w/ater and 10 788 fecal samples). The number of samples collected by 
bird species is as follows: chickens, 11 513; ducks, 1616; quail, 1460; pet birds, 2350; lesser w/histling ducks, 2196; and unidentified wild waterfowl, 725. 
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RESULTS 

Influenza A viruses in Bangladeshi live poultry markets 

From November 2008 to June 2013, 975 low pathogenic avian influ- 
enza (LP AT) H9N2 viruses and 66 1-IPAI 1-I5N1 viruses were isolated 
from 19 897 biological and environmental samples collected from 
apparently healthy birds in Bangladesh. Co-infection with both 
HPAl H5N1 and LPAl H9N2 viruses was found in 78 samples from 
single hosts (Figure 1 ). Influenza viruses from subtypes other than H5 
and H9 were isolated very rarely, mainly H3N8 in one of the duck 
farms in 2012. 

HPAI H5N1 viruses were isolated only from the live bird markets 
from three bird species: chickens, quail and ducks. AIVs were not 
detected in samples from migratory wild waterfowl collected at the 
wetland and the lake. 

During the study period, H9N2 viruses were isolated every month 
except September 2009 and there was no specific pattern of seasonality 
in virus prevalence. H5N1 viruses were mainly isolated during the 
winter and spring seasons of each year, but sporadically one to five 
H5N1 viruses or a co-infection with H5N1 and H9N2 in a single-host 



sample were detected from August to October. Figure 1 shows the 
number of samples collected and viruses isolated every month during 
the surveillance period. 

Antigenic analysis in HI assays 

Table 2 presents results of the HI assays with 30 H5N1 isolates from 
Bangladesh, using the World Health Organization reference antisera 
against H5N1 vaccine candidate strains from clades 2.3.2.1 and 2.2 
and ferret sera produced against the H5N1 viruses isolated from 
Bangladesh. 

The two H5N1 isolates from clade 2.2.2 (A/Ck/BD/9636/10 (H5N1) 
and A/Ck/BD/9675/11 (H5N1)) cross-reacted with antisera to clades 
2.2.2 and 2.3.2.1 viruses in the HI assays. Their HI titers against a-A/ 
Bar-headed goose/Qinghai/lA/05 (H5N1), clade 2.2, were within 1 
log2 difference from the homologous titer of that serum, and were 
as high as the HI titers that the clade 2.2.2 isolates showed against 
the sera produced to the recent Bangladeshi clade 2.3.2.1 isolates ot-A/ 
Ck/BD/15205/12 (H5N1) and ot-A/Env/BD/15121/12 (H5N1). 
In contrast, the ot-A/Bar-headed goose/Qinghai/lA/05 (H5N1), 



Table 2 Antigenic analysis of H5N1 influenza viruses from Bangladesh in the hemagglutination inhibition assays 



H5N1 antigen a-HSNl (post-infection terret antisera) 







a-BHG/ 




a-BS/HK/ 


d-CM/HK/ 


a-Ck/BD/ 


a-Env/BD/ 


a-Dk/BD/ 


Reference antigens 


Clade 


QH/IA 


o!-Hubei/l 


1161* 


5052 


15205' 


15121*'' 


19097*'' 




2 2 


640 


80 


320 


80 


80 


640 


160 


nc A/HI )RFi/i /9ni n 


2 3 2 1 


40 


640 


2550 


320 


640 


5120 


1280 


lAO-AV Drt r\ IN o vVrtLLU Vv/ nUI N L] rXUINLi/ i ID i/iiUiU 


2 3 2 1 


40 


160 


5120 


640 


640 


5120 


1280 


RG-A/COMMON MAGPIE/HONG KONG/5052/2007 


2.3.2.1 


40 


160 


1280 


640 


320 


2560 


640 


A/CHICKEN/BANGLADESH/15205/2012 


2.3.2.1 


40 


80 


1280 


320 


640 


2560 


2560 


A/ENVIRONIVIENT/BANGLADESH/15121/2012 


2.3.2.1 


40 


80 


1280 


320 


320 


2560 


1280 


A/DUCK/BANGLADESH/19097/2013 


2.3.2.1 


40 


160 


2560 


640 


1280 


5120 


2560 


TEST ANTIGENS 


















A/CHICKEN/BANGLADESH/9636/2010 


2.2.2 


640 


160 


320 


80 


640 


640 


320 


A/CHICKEN/BANGLADESH/9675/2011 


2.2.2 


320 


80 


80 


40 


320 


320 


160 


A/CHICKEN/BANGLADESH/12197/2011 


2.3.2.1 


40 


80 


1280 


320 


320 


2560 


640 


A/CHICKEN/BANGLADESH/15078/2012 


2.3.2.1 


40 


80 


1280 


160 


320 


1280 


640 


A/CHICKEN/BANGLADESH/15079/2012 


2.3.2.1 


20 


160 


640 


320 


320 


2560 


640 


A/CHICKEN/BANGLADESH/15083/2012 


2.3.2.1 


40 


160 


1280 


320 


320 


2560 


640 


A/CHICKEN/BANGLADESH/15085/2012 


2.3.2.1 


10 


160 


640 


320 


320 


1280 


640 


A/ENVIRONMENT/BANGLADESH/15105/2012 


2.3.2.1 


10 


40 


1280 


320 


320 


2560 


640 


A/ENVIRONMENT/BANGLADESH/15114/2012 


2.3.2.1 


<10 


80 


1280 


160 


320 


2560 


640 


A/ENVIRONMENT/BANGLADESH/15126/2012 


2.3.2.1 


20 


80 


1280 


160 


320 


5120 


640 


A/ENVIRONIVIENT/BANGLADESH/15145/2012 


2.3.2.1 


10 


160 


2560 


320 


640 


5120 


1280 


A/ENVIRONMENT/BANGLADESH/15147/2012 


2.3.2.1 


20 


80 


1280 


320 


320 


5120 


1280 


A/ENVIRONMENT/BANGLADESH/15165/2012 


2.3.2.1 


10 


160 


5120 


640 


640 


5120 


1280 


A/CHICKEN/BANGLADESH/15210/2012 


2.3.2.1 


20 


160 


640 


80 


640 


2560 


640 


A/ENVIRONIVIENT/BANGLADESH/15661/2012 


2.3.2.1 


10 


80 


1280 


320 


640 


2560 


640 


A/ENVIRONIVIENT/BANGLADESH/16298/2012 


2.3.2.1 


10 


160 


1280 


320 


640 


2560 


1280 


A/DUCK/BANGLADESH/17753/2012 


2.3.2.1 


40 


160 


1280 


320 


320 


2560 


640 


A/CHICKEN/BANGLADESH/18061/2012 


2.3.2.1 


40 


40 


640 


80 


320 


2560 


640 


A/CHICKEN/BANGLADESH/18065/2012 


2.3.2.1 


40 


80 


640 


80 


320 


2560 


640 


A/CHICKEN/BANGLADESH/18248/2012 


2.3.2.1 


10 


160 


1280 


320 


640 


2560 


640 


A/ENVIRONMENT/BANGLADESH/18948/2013 


2.3.2.1 


20 


80 


1280 


320 


640 


2560 


640 


A/DUCK/BANGLADESH/18949/2013 


2.3.2.1 


<10 


160 


2560 


320 


640 


5120 


1280 


A/CHICKEN/BANGLADESH/18984/2013 


2.3.2.1 


20 


40 


640 


40 


320 


2560 


640 


A/ENVIRONIVIENT/BANGLADESH/18989/2013 


2.3.2.1 


20 


40 


320 


40 


160 


1280 


320 


A/DUCK/BANGLADESH/19099/2013 


2.3.2.1 


<10 


40 


640 


160 


320 


1280 


320 


A/QUAII_/BANGLADESH/19250/2013 


2.3.2.1 


<10 


20 


160 


40 


40 


320 


160 


A/QUAII_/BANGLADESH/19254/2013 


2.3.2.1 


40 


10 


160 


40 


40 


320 


80 


A/ENVIRONMENT/BANGLADESH/19338/2013 


2.3.2.1 


<10 


20 


320 


40 


80 


640 


160 



'During serum production, ferrets were boosted witti antigen in incomplete Freund's adjuvant on day 14 post-inoculation. 
' indicates sera that have been produced against recent highly pathogenic H5N1 viruses from Bangladesh. 
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Table 3 Amino-acid residues characteristic of the H5N1 viruses isolated from Bangladesh 





Amino-acid 






H5N1 BD Clade 


H5N1 BD 




Gene segment 


residue 


Specific function 


2.2.2 


Clade 2.3.2.1 


Comments 


HA 


140 


Anti 


genie site 


R 


N 






141 


Anti 


^enic site 


S 


S 






154 


Anti 


^enic site 


N 


D 






155 


Anti 


genie site 


D 


N* 


Increases virus binding to a-2,6 sialic acid receptors and reduces lethality and 




156 


Anti 


genie site 


A* 


A*/T 


systemic spread in mice.^^ A/Ck/BD/15210/12 (2.3.2.1) tiave HA N155D 




162 


Anti 


genie site 


1 


K 


A/Ck/BD/12197/11 (2.3.2.1) has K162R 




163 


Anti 


genie site 


S 


S 


A/Ck/BD/15079/12 and A/Ck/BD/15083/12 have S163D 




189 


Anti 


genie site 


K 


R*/G 


Increases virus binding to o(-2,6 sialic acid receptors.^" A/Ck/BD/15210/12 
(2.3.2.1) has HA R189K 




94 






N* 


N*/T 


Enhances receptor specificity for a-2,6 sialic acid receptors and enhances viral 




227 (226) 






M 


1* 


infectivity of AA/ietnam/1203/04 in mammalian cells^^ 




14 

88 






E 
D 


K* (also N, M) 
G* 


Also in clade 2.3.2.1 isolates from India and Nepal and some Cambodian isolates 

from clade 1.1 
Also in clade 2.3.2.1 isolates from India and Nepal 




115 






Q 


R* 


Also in clade 2.3.2.1 isolates from India and Nepal 




446 






K 


R* 


Also in clade 2.3.2.1 isolates from India and Nepal 


NA 


166 
374 






A 
1 


V* 
V* 


Also in clade 2.3.2.1 isolates from India 
Also in clade 2.3.2.1 isolates from India 


PB2 


102 






N 


S* 


Also in clade 2.3.2.1 isolates from India 
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clade 2.2, antiserum cross-reacted with clade 2.3.2.1 viruses with very 
low HI titers (sC40). 

Clade 2.3.2.1 viruses from Bangladesh reacted against aU clade 2.3.2.1 
antisera. The cross-reactivity was lowest to ot-A/Hubei/1/10 (H5N1) 
(clade 2.3.2.1), with HI titers between 10 and 160, whereas the homo- 
logous titer of the antiserum was 640. The clade 2.3.2.1 viruses had 
highest HI titers against a-A/Env/BD/15121/12 (H5N1) (clade 2.3.2.1) 
and a-A/Barn swaUow/HK/1 161/10 (H5N1) (clade 2.3.2.1), followed by 
a-A/Dk/BD/19097/13 (H5N1) (clade 2.3.2.1). The high HI titers 
obtained with these sera (1/2560-1/5120) were because the immune 
response was boosted with an adjuvant, but this would likely not affect 
the antigenicity pattern of the viruses. The HI titers of the clade 2.3.2.1 
isolates against a-A/Ck/BD/15205/12 (H5N1) (clade 2.3.2.1) were 
within 1 log2 difference of its homologous titer. This observation con- 
firmed that sera specifically raised to recent Bangladeshi isolates pro- 
vided the highest level of HI. Clade 2.3.2.1 viruses from Bangladesh also 
reacted well to ot-A/Common magpie/HK/5052/07 (H5N1) (clade 



2.3.2.1), but the HI titers were up to 4 logS2 lower than the homologous 
titer. Three recent isolates from clade 2.3.2.1— A/Quail/BD/19250/13 
(H5N1), A/Quail/BD/19254/13 (H5NI) and A/Env/BD/19338/13 
(H5N1) — had very low HI titers ( <40 with aU non-boosted sera and 
160-640 with boosted sera), supporting the occurrence of antigenic drift 
among recent H5N1 isolates. 

Molecular characteristics of Bangladeshi H5N1 viruses 

Alignment of the sequences of the H5N1 viruses isolated in this study 
with representative sequences from NCBI and Global Initiative on 
Sharing AU Influenza Data identified some amino-acid residues in each 
genomic segment that were molecular signatures of the viruses isolated 
in the South Asian geographic region encompassing Bangladesh- 
India-Nepal (Table 3). This molecular signature was clade-specific. 

HA gene. AU H5N1 viruses isolated in Bangladesh had multiple basic 
amino acids at the cleavage site (CS) at positions 321-333 of the HA 



Emerging Microbes and Infections 



H5N1 viruses in Bangladesli 

Marinova-Petkova etal 



99i 



99 



99 



99 



86 
98 



A Quail Bangladesh 1 9250 1 3 
A Quail Bangladesh 1 9254 1 3 

A Environment Bangladesh 18948 13 
A chicken Bangladesh 12VIR 7140 11 2012 
r bA Duck Bangladesh 1 7753 1 2 
I — "A Duck Bangladesh 19097 13 

bA Duck Bangladesh 19099 13 

■A Environment Bangladesh 15661 12 
■A Environment Bangladesh 16298 12 
jrA chicken Bangladesh 11rs1984 16 2011 
qlAcrow Bangladesh 11rs1984 11 2011 
■A Chicken Bangladesh 12197 11 

"A Chicken Bangladesh 15078 12 

A duck India 02CA1 0 2011 
_9oA chicken Nepal T9 BH 2011 

I— A chicken Nepal PT 16 2012 
—A chicken Bangladesh 12VIR 7140 1 2011 
A Chicken Bangladesh 15085 12 
ADuckBangladesh1894913 
A chicken India 09CA10 2011 

■A Environment Bangladesh 15147 12 
r»A Environment Bangladesh 15105 12 
IbA" 



~99"1bA Environment Bangladesh 15114 12 
■A Environment Bangladesh 15126 12 
iA Environment Bangladesh 15121 12 
■A Environment Bangladesh 15145 12 
lA Environment Bangladesh 15165 12 
-A chicken India 03CL488 2011 

BAChicken Bangladesh 15210 12 

lA Chicken Bangladesh 18248 12 



93 



98 



98 



A chicken Bangladesh 12VIR 7140 10 2012 
— bA Environment Bangladesh 19338 13 
99|BA Chicken Bangladesh 18061 12 
_J BAChicken Bangladesh 18065 12 
"99[p BAChicken Bangladesh 18984 13 
98LbA Environment Bangladesh 18989 13 
AHubei 1 2010 

duck Hunan S4030 2011 
A chicken Vietnam NCVD 1963 2012 
-A chicken Jiangsu k0402 2010 
99rA grebe Tyva 2 2010 
7<i ^-A whooper swan Mongolia 1 2010 

W A chicken Miyazaki 11 2011 

La great crested grebe Qinghai 1 2009 
A Hong Kong 6841 2010 

-A duck Vietnam NCVD 1547 2012 



J AH 

4Li Ac 

72 A 



Hubei-like 




2.3.2.1 



99 



-Aduck Hunan S4150 2011 HK- 
»A Chicken Bangladesh 15079 12 like 
991bA Chicken Bangladesh 15083 12 
-Aduck Hunan S4220 2011 I 
-A barn swallow Hong Kong 1161 2010 
-Aduck Vietnam NCVD 1463 2012 
-A duck Hunan S4234 2011 hs. 
-A wild duck Fujian 1 2011 [j^g 
-A chicken Lao LH1 2010 '2 34 1 




~i ^A duck Vietnam NCVD 422 201 0 

A chicken Laos P01 30 2007 12.3.4 
A village chicken Kyaing Tong 2433 2007 
Aduck Lao 961 2010 

A chicken Bangladesh 11rs1984 30 2011 
Achicken Bangladesh 11rs1984 37 2011 
99, Achicken Menofia CAI10 2010 




2.3.4.2 



Aduck El Fayoum CAM 2010 

-A duck Egypt Q4596D 2012 



2.2.1.1 



91 



gaA Bar headed Goose Qinghai 59 2005 

A bar headed goose Qinghai 051 0 2005 I 
—A chicken India WB NIV92456 2009 



I 2.2.1 



2.2 



99 



A duck Tripura 1 03597 2008 
— Achicken Bangladesh 1151 11 2010 

'IP A chicken West Bengal 155505 2009 

rn r iA chicken Bangladesh 9636 1 0 
' — rA chicken Bangladesh 11rs1984 4 2010 
A Chicken Bangladesh 9675 11 
A/Bangladesh/3233/2011 
99pA chicken Cambodia 083LC 2011 
^A heron Cambodia TM068 201 1 



9ilr 



2.2.2 



1.1 



97'- 



-A Cambodia P0322095 2005 
A Vietnam 1203 2004 
A duck Hunan 795 2002 12.1.3 

A Indonesia CDC1032 2007 
-Achicken East Java UT6045 2007 
87l — Achicken East Java UT6020 2006 
-A goose Guangdong 1 1996 



1 



I 2.1.3.2 
I 2.1.3.3 
I 2.1.3.1 



I 0 



10 



Figure 2 Phylogenetic tree of tlie HA genes of HPAI H5N1 viruses from Bangladesh, generated by neighbor-joining method In MEGA 5. Numbers at the branches 
Indicate bootstrap values; only values >70 are shown. ▲ Indicates Bangladesh Isolates clade 2.2.2; ■ Indicates Bangladesh Isolates from the A/Hubel/l/10-lll<e 
lineage of clade 2.3.2.1; • Indicates Bangladesh Isolates from the A/Hong Kong/6841/lO-llke lineage of clade 2.3.2.1. 
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(H5 numbering), which is a marker for high pathogenicity in chick- 
ens.^* Clade 2.2.2 isolates had a PQGERRRKKR/GLF CS sequence, 
which was identical with the HA CS motif of wild-bird H5N1 isolates 
from Qinghai Lake in 2005^^ and clade 2.2.2 isolates from the India 
and Bangladesh geographic region, including the human isolate A/ 
BD/3233/2011 (H5N1) (GenBank accession number AEA50985). 

Of the 28 clade 2.3.2.1 isolates, 26 had the PQRERRRKR/GLF CS 
motif, including deletion at position 329, which is common in clade 
2.3.2.1 H5N1 viruses from the Bangladesh, India and Nepal regions 
and from wild birds from Qinghai (China), Mongolia, and Tyva in 
2009 and 2010. Two isolates from clade 2.3.2.1 had an amino-acid 
variation at position 323 of the CS: A/Ck/BD/15210/12 (H5N1) had 
R323G and A/Env/BD/19338/13 (H5N1) had R323K. 

Table 3 lists the amino-acid residue positions at the HA antigenic 
sites of isolates from clades 2.2.2 and 2.3.2.1. Residue R189G found in 
several clade 2.3.2.1 isolates from Bangladesh was unique, whereas the 
others residues were common in H5N1 viruses from Bangladesh, India 
and Nepal. HA residues 94N, 155N, 156A, 189R and 2261, which we 
found to occur in the majority of H5N1 viruses isolated in Bangladesh, 
have been reported to increase virus binding to a-2,6 sialic acid recep- 
tors (Table 3).^^^"^^ Residues 222Q and 224G of aU Bangladeshi H5N1 
isolates remained unchanged, suggesting an a-2,3-linked sialic acid 
preference and a greater likelihood of avian infectivity.^* In the HA 
receptor-binding site, isolate A/Ck/BD/15210/12 (H5N1) had a single 
change within the 220-loop (217-224) K218R. 

In the HA of the clade 2.2.2 isolates from Bangladesh, there were six 
glycosylation sites at positions 10 (NNST), 23 (NVT), 165 (NNT), 286 
(NSS), 484 (NGS) and 543 (NGS) (H5 numbering). In addition to 
these glycosylation sites, clade 2.3.2.1 isolates also had another glyco- 
sylation site at position 140 (NSS). Isolate A/Ck/BD/ 18248/ 12 (H5N1) 
had loss of a glycosylation site at position 165, which was caused by the 
amino-acid substitution T167I. 

NA. All the Bangladeshi H5N1 viruses had a 20-amino-acid deletion 
(positions 49-68) in the NA stalk region, which is reported to be 
required for influenza viruses to adapt from wild aquatic birds to 
domestic chickens^' and for enhanced virulence in mice.^" Amino 
acids 116V, 119E, 136Q, 156R, 199D, 2231, 247S, 275H, 277E and 
295N (N2 numbering) in the NA conferred sensitivity to oseltamivir 
and/or zanamivir in the Bangladeshi H5N1 isolates.'''"^* Isolate A/Dk/ 
BD/18949/12 (H5N1) had the amino-acid substitution I117V, which 
may cause reduced susceptibility to oseltamivir.^^ Most of the isolates 
had NA 149V, suggesting sensitivity to zanamivir,^^ but we also found 
two mutations 1491 (A/Ck/BD/ 18248/ 12 (H5N1)) and 149F (A/Env/ 
BD/19338/13 (H5N1)) at that position that had not been previously 
reported to cause resistance. 

PB2. The PB2 E627K mutation, which is important for mammalian 
host adaptation in AIVs,^'' was present in the H5N 1 isolates from clade 
2.2.2, but not in any of the viruses from clade 2.3.2.1 isolated in 
Bangladesh, India and Nepal. The presence of amino-acid substitu- 
tions L89V, G309D, R477G, I495V and A676T of the PB2 protein in all 
H5N1 viruses isolated in this study suggests enhanced polymerase 
activity and increased virulence in mice.^^ 

The PB2 amino-acid residue 702R, which was present in most of the 
H5N1 viruses from Bangladesh and India isolated in 201 1 or later (clade 
2.3.2.1) is a molecular marker of human HlNl, H2N2 and H3N2 and 
classical swine isolates and was previously found in the 1918 pandemic 
virus, A/PR/8/34 (HlNl) and seasonal human isolates.^* Isolates of 
avian and equine origin usually possess PB2 amino-acid residue 702K. 



PBl. The isolate A/Dk/BD/ 18949/ 12 (H5N1) had the PBl amino- 
acid substitution I368V. This mutation may be important for trans- 
mission of H5 viruses in ferrets,^ and has been observed in the PBl 
protein sequence of most of the newly emerged H7N9 viruses in east- 
ern China and H9N2 viruses from China.'^ 

PA. Another human-like mutation found in clade 2.3.2.1 H5N1 
avian isolates from Bangladesh and India was PA S421I, which 
has been previously reported in A/Vietnam/1203/2004 (H5N1). In 
A/Ck/Vietnam/NCVD5/2003 (H5N1), this mutation in combination 
with PB2 E627K causes increased pathogenicity, higher replication 
titers and higher mortality in mice.*" 

M. The Ml proteins of isolates A/Ck/BD/ 15079/ 12 (H5N1) and 
A/Ck/BD/15083/12 (H5N1) possess several hallmark amino-acid resi- 
dues of the recently emerged human and avian H7N9 viruses and of 
H9N2 viruses from China (L46I, T140A, V142G, I219V, A227T, 
K242N and M248L). None of these residues was present in the 
Bangladeshi H9N2 isolates. In the M2 protein, isolate A/Ck/BD/ 
15078/12 (H5N1) had the V27A mutation and isolates A/Ck/BD/ 
15079/12 (H5N1) and A/Ck/BD/15083/12 (H5N1) had the S31N 
mutation. Both mutations are known to cause reduced susceptibility 
to amantadine and rimantadine.*' 

NS. AU Bangladeshi H5N1 isolates fr'om clades 2.2.2 and 2.3.2.1 had a 5- 
amino-acid deletion (positions 80-85) in the NS protein, which is also 
present in other HPAI H5N1 viruses isolated after 2001 and is associated 
with increased virulence of H5N1 viruses in chickens and mice.*^ 

NP. Amino-acid changes in the NP genes of the Bangladeshi H5N1 
viruses from clade 2.3.2.1 found also in isolates from India are shown 
in Table 3. 

Phylogenetic analysis of Bangladeshi H5N1 isolates 

To track the evolutionary relationships of the H5N1 viruses isolated 
from Bangladesh with viruses representing different H5N1 clades and 
lineages, as well as other influenza viral subtypes, clade affiliations of 
aU viruses included in the study were assigned in the HA phylogenetic 
tree and then used in the phylogenetic trees of the other seven genes. 

HPAI H5N1 viruses from clade 2.2.2. Phylogenetic analysis of the 
HA gene of 30 H5N1 isolates from Bangladesh showed that two 
isolates— A/Ck/BD/9636/ 10 (H5N1) and A/Ck/BD/9675/11 (H5N1)— 
belonged to clade 2.2.2 (Figure 2). In the phylogenetic trees of all eight 
genes both isolates clustered together and were closely related to H5N1 
viruses from clade 2.2.2 isolated in Bangladesh in 2010 and 2011, 
which all descended from H5N1 isolates from West Bengal and 
Tripura in India from 2008 and 2009. Their common ancestors were 
H5N1 isolates from bar-headed geese from Qinghai Lake in 2005 
(Figures 3-5 and Supplementary Figures S1-S4). The closest relative 
of isolate A/Ck/BD/9675/11 (H5N1) in all genes was isolate A/ 
Bangladesh/3233/2011 (H5N1), which showed 99% similarity in the 
basic local alignment search tool (BLAST) algorithm analysis*^^ in 
NCBI. 

HPAIH5N1 viruses from clade 2.3.2.1. Of the 30 H5N1 isolates in 
this study, 28 (93%) fell within clade 2.3.2.1 of the HA phylogenetic 
tree and 26 (86.7%) clustered together in one distinct group of the 
phylogenetic tree, which also included viruses from India identified in 
201 1 and from Nepal identified in 201 1 and 2012. This monophyletic 
cluster did not change as a group in the phylogenetic analysis of the 
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Figure 3 Pliylogenetic tree of the NA gene of H5N1 viruses from Bangladesh, generated by neighbor-joining method in MEGA 5. Numbers at the branches indicate 
bootstrap values; only values >70 are shown. ▲ indicates Bangladesh isolates clade 2.2.2; ■ indicates Bangladesh isolates from the A/Hubei/l/lO-like lineage of 
clade 2.3.2.1; # indicates Bangladesh isolates from the A/Hong Kong/6841/lO-like lineage of clade 2.3.2.1. 
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Figure 4 Phylogenetic tree of the PB2 gene of H5N1 viruses from Bangladeshi, generated by neigfibor-joining method in MEGA 5. Numbers at tfie branches indicate 
bootstrap values; only values >70 are shown. ▲ indicates Bangladesh isolates clade 2.2.2; ■ indicates Bangladesh isolates from the A/Hubei/l/lO-like lineage of 
clade 2.3.2.1; • indicates Bangladesh isolates from the A/Hong Kong/6841/lO-like lineage of clade 2.3.2.1. 
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eight gene segments and showed evolutionary distance from the clos- 
est relatives outside the group. The cluster shared a common node 
with the A/Hubei/l/lO-like genetic group in the HA, PBl and NS trees 
(Figure 2 and Supplementary Figures SI and S4). In the PA gene tree 
(Supplementary Figure S2), the Bangladeshi clade 2.3.2.1 shared a 
node with a group of clade 2.3.2.1 viruses, which included all repre- 
sentatives of the A/Hubei/l/lO-like lineage as well viruses of the A/ 
HK/6841/lO-like and the A/Barn swaUow/HK/1161/lO-like hneages. 
In the PB2 tree (Figure 4), the big cluster of Bangladeshi clade 2.3.2.1 
viruses shared a common node with the A/HK/6841/lO-like lineage of 
clade 2.3.2.1. Interestingly, in the NA (Figure 3) and NP gene 
(Supplementary Figure S3) trees, the Bangladeshi clade 2.3.2.1 cluster 
was a descendent of A/Dk/Lao/961/10 (H5N1), clade 2.3.4.2, and was 
clearly distinct from all 3 lineages of clade 2.3.2.1 (A/Hubei/l/lO-like, 
A/HK/6841/lO-like and A/BS/HK/1161/lO-like lineages). 

In the M gene tree (Figure 5), the Bangladeshi cluster of clade 2.3.2. 1 
viruses also included A/Dk/Lao/961/10 (H5N1), clade 2.3.4.2. In the 
phylogenetic tree built using the neighbor-joining method, the cluster 
shared a node with the clade 2.3.2.1 viruses from the A/HK/6841/10- 
like lineage (Figure 5). However, in the tree built using the maximum 
likelihood method (not shown), the cluster shared a node with clade 
2.3.4.2. In both cases, clade 2.3.4 was the ancestor of the Bangladeshi 
H5N1 2.3.2.1 viruses with regard to the M gene. Thus, the resultsofthe 
phylogenetic analysis suggest that the majority of the currently cir- 
culating H5N1 viruses from clade 2.3.2.1 in Bangladesh have under- 
gone intraclade reassortment in the PB2 gene with 2.3.2.1 H5N1 
viruses genetically close to the A/HK/6841/lO-like lineage, and inter- 
clade reassortment in the NA, NP and M genes with H5N1 viruses 
from clade 2.3.4.2. 

Since this cluster of Bangladeshi clade 2.3.2.1 viruses is very homo- 
genous (i.e., it always has the same virus pattern in aU trees) and shows 
that all isolates included have the same ancestors, this may be a sign of 
a single introduction of clade 2.3.2.1 in 2011 into the geographic 
region of Bangladesh, India and Nepal and ongoing evolution within 
the region that is supported by poultry movement. In all the phylo- 
genetic trees in this study, some taxa within the Bangladeshi 2.3.2.1 
cluster were not strongly supported by the bootstrap value (<70%), 
which is indicative of an ongoing evolutionary process. 

Within the HA phylogenetic tree, isolates A/Ck/BD/ 15079/ 12 
(H5N1) and A/Ck/BD/15083/12 (H5N1), both original samples col- 
lected on the same day in lanuary 2012 from the same live bird market 
in Dhaka city, fell within clade 2.3.2. 1, but did not cluster together with 
the rest of the H5N1 isolates from Bangladesh. Both isolates had 
identical HA nucleotide sequences and were genetically related to A/ 
duck/Hunan/S4150/ll (H5N1). They all feU within the A/Hong Kong/ 
6841/2010 (H5Nl)-like genetic lineage of clade 2.3.2.1. This genetic 
group mainly has H5N1 viruses from Japan, China, Vietnam and 
Indonesia (Figure 2). 

Detection of reassortant H5N1 viruses with an H9N2-like M gene 

The A/HK/6841/10 (H5Nl)-like viruses A/Ck/BD/15079/12 (H5N1) 
and A/Ck/BD/15083/12 (H5N1) were closely related to A/duck/ 
Hunan/S4150/ll (H5N1) in the H5, Nl, PB2, PA, NP and NS gene 
phylogenetic trees and closely related to A/duck/Hunan/S4150/ll 
(H5N1) and A/duck/Hunan/S4220/ll (H5N1) in the PBl phyloge- 
netic tree. However, phylogenetic analysis of the M gene revealed that 
the M genes of isolates A/Ck/BD/15079/12 (H5N1) and A/Ck/BD/ 
15083/12 (H5N1) were derived from H9N2 viruses isolated in eastern 
China and are closely related to A/Chicken/Wuxi/7/10 (H9N2) and 
A/Chicken/Zhejiang/329/11 (H9N2) (Figure 5). Nucleotide BLAST 



analysis in NCBI of the M gene sequence of isolates A/Ck/BD/ 15079/ 
12 (H5N1) and A/Ck/BD/15083/12 (H5N1) showed maximum identity 
of 99% with the M gene of A/Chicken/Rizhao/1313/13 (H9N2) and 
some other H9N2 viruses isolated in 2013 from the same location, as 
well as A/Chicken/Hong Kong/JV75/ll (H9N2), A/Chicken/Zhejiang/ 
329/11 (H9N2), and isolates from Shanghai and other locations in east- 
ern China. 

Many H9N2 viruses were isolated from the same Uve bird markets in 
Dhaka city during the last few years, but we did not detect reassort- 
ment between the H5N1 and H9N2 lineages circulating in Bangladesh. 
The Bangladeshi H9N2 isolates represented the Gl genetic lineage and 
in the M gene phylogenetic tree, they fell within an evolutionarily 
distinct cluster from the H9N2 isolates from China (Figure 5). The 
M gene phylogenetic tree also revealed that the M gene of isolates A/ 
Ck/BD/15079/12 (H5N1) and A/Ck/BD/15083/12 (H5N1) and of 
recently emerged H7N9 viruses isolated from humans and birds in 
China have common ancestors — M genes of AI H9N2 viruses circulat- 
ing in China. Negovetich et alP reported in 2011 that influenza 
viruses from subtypes HI, H3, H4 and HIO have been isolated from 
poultry markets in Bangladesh and we detected LPAI H3N8 in a duck 
farm in 2012; however, reassortment between any of the abovemen- 
tioned LPAI subtypes and H5N1 was not detected. 

DISCUSSION 

Bangladesh experienced its first HPAI H5N1 outbreak in February 
2007, and the isolated viruses belonged to clade 2.2 of the Qinghai 
lineage that emerged in 2005.'''' Phylogenetic analyses'''' and evidence 
of wild-bird migration''^ suggested that the virus had been introduced 
into Bangladesh by migratory birds from Mongolia in early 2007. 
Since then, the HPAI H5N1 virus has become endemic to the poultry 
in Bangladesh. 

In January 2011, clade 2.3.2.1 was introduced into the country. 
Islam et al." reported the first cases of the new clade in a crow, a duck 
and a quail, and Khan et al^^ reported a crow die-off caused by clade 
2.3.2.1 in Bangladesh, which had started on 17 January 2011. 

Phylogenetic analyses of the earliest clade 2.3.2.1 isolates from 
Bangladesh revealed a close relationship with H5N1 viruses from 
Nepal from 2010 and India from February 201 1, with wild birds from 
Mongolia and Tyva from 2009 to 2011 being the common ancestors 
for all these viruses. "'''^ Although the exact pathway of introduction of 
clade 2.3.2.1 into Bangladesh remains unclear, it is possible that virus 
transmission to neighboring territories occurred through crows, after 
the virus had been already introduced into Nepal from Mongolia and 
Tyva via migratory birds previously. Clade 2.3.2.1 is currently circulat- 
ing in poultry in Bangladesh and has displaced clade 2.2.2. 

H5N1 viruses from clade 2.3.4.2 were isolated in Bangladesh in 
February and March 2011,'"'"'''^ but these viruses were only sporadi- 
cally detected and did not establish a sustained circulation. 

We report here the isolation of 66 HPAI H5N1 viruses and further 
analyses of 2 H5N1 viruses from clade 2.2.2 from December 2010 and 
January 2011 and of 28 H5N1 isolates from clade 2.3.2.1 from 2011 to 
2013. Until February 201 1, we only detected H5N1 viruses from clade 
2.2.2. During March-July 201 1, we did not detect any H5N1 viruses in 
the samples from Bangladesh, and from August 20 11 to 2013 the only 
H5N1 clade detected was 2.3.2.1, indicating that the clade change 
occurred during the latter period. 

AU HPAI H5N1 viruses and 99% of the H9N2 viruses isolated in this 
study were detected at the live bird markets in Dhaka city, which 
underlines the important role that the markets play in perpetuating 
and maintaining both subtypes. We did not detect H5N1 viruses on 
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Figure 5 Phylogenetic tree of the M gene of H5N1 viruses from Bangladesh, generated by neighbor-joining method in MEGA 5. Numbers at the branches indicate 
bootstrap values; only values >70 are shown. ▲ indicates Bangladesh isolates clade 2.2.2; ■ indicates Bangladesh isolates from the A/Hubei/l/lO-like lineage of 
clade 2.3.2.1; 9 indicates Bangladesh isolates from the A/Hong Kong/6841/lO-like lineage of clade 2.3.2.1. 
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the farms, which may not be representative of all farms in the country. 
During the five years of surveUlance, we did not observe any clinical 
signs of a disease at the farms and retail markets in Bangladesh. 
Normal bird mortality rates were reported by poultry workers and 
retailers at the farms and markets where we were collecting samples. 
The lack of morbidity and mortahty among poultry is unusual for 
HPAI H5N1 viruses, and may be explained with protection against 
the HPAI H5N1 viruses provided by the year-round circulating LPAI 
H9N2 viruses. However, this hypothesis has not been tested because it 
was beyond the scope of our study. 

In January 2012, we detected two reassortant (7:1) HPAI (H5N1) 
clade 2.3.2.1 viruses that were genetically distinct from all the other 
viruses isolated from Bangladesh and possessed the M gene from 
the H9N2 subtype, which circulated in China in 2010-2012. The 
remaining seven genes of isolates A/Ck/BD/15079/12 (H5N1) and 
A/Ck/BD/15083/12 (H5N1) were closely related to A/HK/6841/ 
10-like A/duck/Hunan/S4150/ll (H5N1) fi-om clade 2.3.2.1. The 
phylogenetic analysis revealed that not only the reassortant H5N1 
viruses from Bangladesh but also A/duck/Hunan/S41 11/11 (H9N2) 
and A/duck/Hunan/S4234/ll (H5N1) fall within the Chinese H9N2 
cluster of the M gene tree. All the three above-mentioned duck 
viruses from Hunan were isolated in November 2011 from a back- 
yard-style duck farm with an open environment and no biosecurity 
measures that was located in the Dongting Lake region of China.''^ 
This lake is the third biggest lake in China and is on the East Asian- 
Australasian Flyway of migratory birds. Therefore, on the basis of 
the finding that none of the genes from reassortant H5N1 viruses 
originated from local isolates, we hypothesize that these viruses had 
been directly introduced into Bangladesh either via migratory birds, 
using overlapping flyways from the Dongting region in China to 
Bangladesh, or via poultry trade and movement, combined with 
infected wild-bird migration. This reassortant 7:1 H5N1 genotype 
was detected only once in Bangladesh, representing transient infection 
in chickens. Our surveillance results provide further confirmation that 
the H9N2 viruses circulating in China have been involved in multiple 
reassortment events with both LPAI and HPAI viruses, but with dif- 
ferent outcomes. Although the reassortant viruses reported in this 
study did not affect the H5N1 situation in Bangladesh, the newly 
emerged H7N9 viruses in China caused many human deaths in 2013. 

Another important finding of this study is that despite many years 
of cocirculation of H5N1 and H9N2 viruses in Bangladesh, we did not 
find reassortant H5N1 viruses that possessed internal genes from the 
local H9N2 lineage. This finding is surprising, because in this study 
H5N1 viruses were more often detected as a co-infection with H9N2 
than alone (66 H5N1 isolates versus 78 isolates containing both H5N1 
and H9N2 in a single-host sample). Shanmuganatham et al}^ also 
found no H9N2 viruses from Bangladesh that possessed internal genes 
from H5N1. As of November 2013, one study has reported two reas- 
sortant H5N1 viruses with the H9N2 PBl gene in poultry in 
Bangladesh.''^ Thus, reassortant viruses between the H5N1 and 
H9N2 subtypes apparently occur in Bangladesh, but our surveillance 
suggests that they do not have the fitness to supplant the dominant 
clade 2.3.2.1. This may be because the H9N2 lineage that currently 
circulates in Bangladesh is a reassortant virus possessing three gene 
segments from subtype H7N3.'* 

Our molecular and phylogenetic analyses of the H5N1 viruses 
showed that the majority of H5N 1 isolates from Bangladesh from both 
clade 2.2.2 and clade 2.3.2.1 were genetically similar to and clustered 
together with contemporary H5N1 isolates from India and Nepal. 
The viruses from clade 2.2.2 did not show any evidence of genetic 



reassortment. They seemed to have established a stable lineage for 
the period of their circulation in Bangladesh (February 2007- 
February 2011) without reassorting with other subtypes. 

Of the H5N1 viruses from clade 2.3.2.1, 93% had common ancestry 
and formed an evolutionary distinct 'Bangladesh' cluster in all eight 
gene phylogenetic trees. By conducting full genome phylogenetic ana- 
lyses, we could detect intra- and inter-clade reassortment events that 
occurred in the Bangladeshi clade 2.3.2. 1 isolates and those from India 
and Nepal. The HA, PBl, NS and PA genes had common ancestors 
with the A/Hubei/l/lO-like (H5N1) 2.3.2.1 lineage, whereas the PB2 
genes were similar to A/HK/6841/lO-like (H5N1) 2.3.2.1 viruses. Even 
more interesting was the finding that the NA, NP and M genes were 
clearly related to clade 2.3.4.2. This clade was detected in the 
Chittagong district in East Bangladesh in February-March 201 1 when 
clade 2.3.2.1 had already been introduced into the country, but it 
circulated for a short period of time and did not spread to other 
districts. Because there is limited sequence information for clade 
2.3.4.2 from Bangladesh in public databases, we could not conclude 
whether clade 2.3.4.2 from Bangladesh had been involved in the reas- 
sortment events of clade 2.3.2.1, especially with regard to internal 
genes. In the NA tree, the only clade 2.3.4.2 virus from Bangladesh 
(A/Ck/Bangladesh/BL/543/11 (H5N1)) seems to be genetically related 
to viruses from Kyang Tong (Myanmar), Vietnam and Laos, but in a 
cluster that is distinct from A/Dk/Lao/961/10, which shares a node 
with the Bangladeshi clade 2.3.2.1 isolates. Thus, it is possible that the 
H5N1 clade 2.3.2.1 had been introduced in the region of Bangladesh, 
India and Nepal after the reassortment events involving clade 2.3.4.2 
occurred. This theory would explain why the majority of isolates from 
the region have the same gene constellation. 

Our overall results from the phylogenetic analyses suggest that both 
clades 2.2.2 and 2.3.2.1 might have had one initial introduction each 
into the region by migratory birds, followed by an intensive evolu- 
tionary process driven by poultry movements and cross-border viral 
transmission among Bangladesh, India and Nepal. New H5N1 viruses 
were introduced into Bangladesh either via wild-bird migrations or by 
poultry trade, but they have not persisted in the poultry population. 
The molecular analysis of the H5N1 viruses isolated in Bangladesh 
showed that they share many common molecular changes in all eight 
genes with isolates from India and Nepal, forming a clade-specific 
genetic signature of the H5N1 viruses that evolved within this region. 
This is valid for both clades 2.2.2 and 2.3.2.1. 

Also, some human-like mutations were found in the genome of the 
H5N1 viruses circulating in Bangladesh, which might explain the low 
number of human H5N1 cases reported in the country. In contrast to 
clade 2.2.2, the currently circulating clade 2.3.2.1 does not carry the PB2 
E627K mutation, which is important for mammalian host adaptation. 
The pathogenicity role of the two human-like mutations PB2 K702R and 
PA S421I, which are a part of the genetic signature of the clade 2.3.2.1 
H5N1 viruses from Bangladesh and surrounding regions and have pre- 
viously been described in the 1918 pandemic virus and in the AA'^ietnam/ 
1203/2004 (H5N1), respectively, needs to be further investigated. 

Several mutations found in the Bangladeshi isolates have been pre- 
viously reported to cause binding to a-2,6 receptors and increased 
pathogenicity in mice.^'"^^ Molecular analysis showed that the NA 
protein of one of the recent isolates had the II 17V mutation, which 
may cause reduced susceptibility to the commonly used neuramini- 
dase inhibitor oseltamivir. Also, three other viruses from 2012, includ- 
ing the reassortant H5N1 isolates, had M2 protein mutations V27A 
and S3 IN, which can cause reduced susceptibility to the M2 ion chan- 
nel blockers amantadine and rimantadine. 



Emerging IVIicrobes and Infections 



H5N1 viruses in Bangladesh 

Marinova-Petkova etal 



Our HI assays showed that all tested H5N1 isolates from the currentiy 
circulating Bangladeshi clade 2.3.2.1 isolated until 2013 are antigenically 
homogenous. Further surveillance activities are required to clarify 
whether the lower titer of the most recent clade 2.3.2.1 isolates is due 
to antigenic drift. The variety of host bird species (chickens, ducks and 
quail) of the H5N1 viruses in Bangladesh did not lead to antigenic 
diversity. Although the HA glycoprotein of the majority of the currently 
circulating clade 2.3.2.1 viruses is genetically closely related to that of the 
Hubei/l/lO-like lineage of clade 2.3.2.1, antigenically, they do not show 
cross-reactivity with this virus. Because of the poor reactivity of the 
Bangladeshi Hubei/l/lO-like clade 2.3.2.1 viruses to the prototype virus, 
the World Health Organization recentiy chose the isolate A/duck/ 
Bangladesh/ 19097/20 13 (H5N1) reported in this study as a new candid- 
ate vaccine virus to represent the H5N1 lineage circulating in this region.^ 

The results of our study emphasize the need of continuous surveil- 
lance activities in Bangladesh that would allow detection of any newly 
emerged AlV with pandemic potential. 
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